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ABSTRACT 

Viscous fluid flows with curved streamlines can support both centrifugal and viscous 
travelling wave instabilities. Here the interaction of these instabilities in the context of the 
fully developed flow in a curved channel is discussed. The viscous (Tollmien-Schlichting) 
instability is described asymptotically at high Reynolds numbers and it is found that it 
can induce a Taylor-Gortler flow even at extremely small amplitudes. In this interaction, 
the Tollmien-Schlichting wave can drive a vortex state with wavelength either comparable 
with the channel width or the wavelength of lower branch viscous modes. The nonlinear 
equations which describe these interactions are solved for nonlinear equilibrium states. 


This work was supported by the National Aeronautics and Space Administration under NASA Contract 
Nos. NAS1-18107 and NASl-18605 while the second and third authors were in residence at the Institute for 
Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, 
VA 23665. 


1. Introduction. 

There are many fluid flows of practical importance where destabilizing centrifugal and 
viscous instability mechanisms are both present. Thus, for example, the flow in a curved 
rectangular duct or the flow over parts of a laminar flow wing (see Harvey and Pride 
(1982)) can support both Taylor-Gortler vortices and Tollmien-Schlichting waves. Indeed 
the latter flow can in turn support Rayleigh instabilities associated with the highly in- 
flectional velocity profiles in some directions. In these flows the possibility exists that 
the nonlinear interaction of the different instability mechanisms might produce premature 
transition to turbulence. In that regard a potentially significant result from the recent 
work of Hall and Smith (1988a) is that interacting oblique Tollmien-Schlichting waves can 
generate longitudinal vortices essentially identical to Taylor-Gortler vortices even in the 
absence of wall curvature. 

One of the first calculations on the interaction of Tollmien-Schlichting (TS) waves and 
Gortler vortices was given by Nayfeh (1981). Nayfeh discussed the effect of a Gortler 
vortex of given size on the growth of oblique TS waves. The amplitude of the vortex was 
assigned arbitrarily however and an eigenfunction shape found by solving the parallel-flow- 
approximation Gortler-vortex equations. It is now known (Hall 1982a, 1983) that, in the 
wavenumber regimes considered by Nayfeh, these equations have solutions of no relevance 
to spatially growing vortices. Furthermore, the amplitude of the vortex cannot be assigned 
arbitrarily; it must of course be determined by either a numerical or analytical solution of 
the Navier-Stokes equation as in Hall (1982b, 1988). Having made these approximations 
Nayfeh found that the vortices could have a massive effect on the growth of TS waves. 
Later Malik (1988) showed that Nayfeh’s numerical calculations were incorrect and that 
his conclusions were in error. 

More recently Bennett and Hall (1988) examined the effect of finite amplitude Gortler 
vortices in fully developed flows on the growth of lower branch TS waves. Here the non- 
linear vortex state was found by solving the Navier-Stokes equations and a linear stability 
analysis of the solutions was given. It was shown that even small amplitude vortices can 
have a significant effect on TS growth rates. The vortex flow found by Bennett and Hall 
had a wavelength comparable to the depth of the channel in which the flow occurred. The 
asymptotic structure of TS waves appropriate to the lower branch of the neutral curve 
corresponds in general to TS waves with a small spanwise wavenumber and Bennett and 
Hall showed how this structure could be modified to allow for a faster spanwise dependence 
induced by the vortices . 

The first nonlinear description of the interaction of vortices and TS waves was given 
by Hall and Smith (1988a). Here it was shown that long wavelength vortices and oblique 
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TS waves undergo a resonant triad interaction at small amplitudes. The interaction is, 
in the first instance, governed by ordinary differential triad amplitude equations which 
possess a finite time singularity. At higher amplitude the interaction is controlled by 
a coupled partial differential-ordinary integro-differential equation system. The solution 
of the system was found to depend crucially on the orientation of the TS waves to the 
vortices. The most dangerous type of interaction concerns TS waves propagating in a 
direction making an angle of more than 41.6° to the vortices since the resulting interaction 
produces a finite time singularity. 

In this paper we shall investigate the stronger type of interaction which occurs after 
this singularity develops. The size of the TS waves is fixed by the condition that the 
downstream velocity component is perturbed at zeroth order by the vortex induced by 
the interacting TS waves. Thus this interaction leads to an 0(1) change from the basic 
state that exists in the absence of vortices. If the size of the TS wave is decreased then 
a weakly nonlinear bifurcation governed by a cubic order amplitude equation typical of 
those found when using the Stuart-Watson method is retrieved. Surprisingly the stronger 
type of interaction can occur at extremely small TS amplitudes and both with and without 
wall curvature being present. The fact that such a relatively small TS wave can have such 
a large effect on the basic state is due to the large initial forces associated with the small 
streamwise length scale of the waves. The nonlinear equilibrium states appropriate to 
this interaction must be found numerically. A similar extension of the Bennett and Hall 
work into the nonlinear regime is also discussed for vortex flows of cross-stream wavelength 
comparable to the channel depth. We note here that other nonlinear states are accessible 
by routes other than that described in this paper; see, for example, Hall and Smith (1988b) 
where the vortex-TS interaction for external flows is discussed. The question of whether 
the present route described is physically the most relevant or whether one of the ’by-pass’ 
routes dominates in an experiment cannot yet be answered. 

In addition, at small vortex wavenumbers the three-dimensional breakdown of a 2D-TS 
wave can be described by an analysis of our interaction equations. Thus we determine the 
size of 2D-TS waves which are neutrally stable to oblique TS waves. A key feature of this 
secondary instability process is the longitudinal vortex system induced by the interacting 
oblique waves. Furthermore the interaction can occur in a straight channel, thus yielding 
a mechanism for the breakdown of 2D-TS waves in parallel flows. 

The procedure adopted in the rest of this paper is as follows: in §2 the nonlinear 
equations governing fully nonlinear vortex flows in curved channels are described; in §3 
these equations are solved in the presence of vortices of 0(1) cross-stream wavenumber 
and the possible equilibrium states are described; in §4 the corresponding calculation for 
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small wavenumber vortices is described; and in §5 the instability of 2D-TS waves in a 
straight channel is discussed through the interaction equations derived in §4. Finally in §6 
we discuss our results and draw some conclusions, mainly that nonlinear effects lead to a 
supercritical bifurcation to a mixed vortex-TS state. 

2. Formulation of the disturbance equations for curved channel flows. 

Consider the flow of a viscous incompressible fluid in a curved channel with walls defined 
by r* = a,a + d with respect to cylindrical polar co-ordinates (r*, 0*, z*). It is assumed 
that the curvature parameter, 6, defined by 

6 = -, ( 2 . 1 ) 

a 

is small. The flow is driven by the streamwise pressure gradient 

l dp* _ -12uV m r* , 00 , 

pd0* d? 1 } 

where p and v are the density and kinematic viscosity of the fluid respectively, whilst V m 
is a typical streamwise flow speed. The pressure gradient (2.2) drives a velocity field v m 
in the streamwise direction and 

v m = V m U{y), (2.3) 


where for small values of 6 


U(y) = U 0 (y) + 0(«), 


U 0 (y) = 6y(l - y), y 


(2.5 a, b) 


The dimensionless variables x, z, and t are then defined by 

aO* z* t*u , , x 

x = — , ,z=—, t = — . (2.6a, 6 , c) 

where t* denotes the (dimensional) time. If (u, v, w) is the velocity field scaled on V m with 
respect to (x,y, z) and the pressure is scaled on pV^ then the Navier-Stokes equations can 

be written in the form 

1 du 6u dv dw _ 

J dx "I" / dy **" dt ’ 


1 ry-,1 * fi 2 1 26 dv 1 dp 6uv 

^{V - a ,--}u+— 2rx - J - = Nu+—, 

1 Cr~r2 a ^ 1 2(5 du d P _ at ^ 

dt J 2 > v RePdx dy~ V 7’ 


(2.7a, b, c,d) 
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where the Reynolds number Re is defined by 


Re = V m - 


( 2 . 8 ) 


and 


J = 1 + By, 


V 2 =jidl + dl + jd v + dl 

N = (u/7)d x + vd y + wd g . 
The Taylor number T is then defined by 

T = 4 Re 2 6, 


(2.86, c,d) 


(2.9) 


and it is known from the work of Dean (1928) that instability in the form of Taylor vortices 
occurs first for 0(1) spanwise wavelengths with T also an 0(1) quantity. In view of the 
smallness of 6 it follows from (2.9) that Taylor vortices first occur at high values of the 
Reynolds number. The flow is therefore susceptible to lower branch Tollmien-Schlichting 
instabilities with wavenumbers of order Re~ j in the azimuthal direction. It is the nonlinear 
interaction of these two modes of instability which will be discussed in the next section. 

However, because a primary aim of the present calculation is to shed light on the related 
external flow problem where long wavelength Gortler vortices are important in the initial 
development of the flow, we consider a closely related interaction problem arising from the 
weakly nonlinear theory of Hall and Smith (1988a). Here the Taylor vortices occur with 
axial wavenumbers 0 [Re~ y) which are comparable with the wavenumbers of lower branch 
Tollmien-Schlichting waves. These long wavelength vortices occur at relatively high values 
of the Taylor number T ~ 0 (Ref) which requires that Re ~ 0(6“ &) whereas the previous 
interaction has Re ~ 0(6“ *). Thus, for a given value of Re, the former type of interaction 
occurs first when 6 is increased from zero. Alternatively, for a given channel with 6 fixed, 
the former interaction occurs at the lowest values of Re. However there is some indication 
from the work of Hall and Smith (1988a, b) that the second interaction can occur even at 
zero values of T, and so it is important to understand its structure. 

Finally in this section we write down the form of the solutions of (2.7) appropriate 
to fully nonlinear Taylor vortices. Firstly, for 0(1) vortex wavelengths the pressure p and 
velocity field (u, v,w) can be written as 


P = 


12x 


+ 


Re 2 Re 2 


][l + 0(Re )] 


(2.10a, 6) 


(u,v,w) = [U 0 + u,-^-,^-][l + 0{ R e 
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Here u,v,w,p are functions of t, y, z and satisfy 


dv dw 

a? + al = 0 

[V\ - d t ]u - \vU' Q = Uvu y + wv z ) 


[^2 - dt\v - ^ ^ uU 0 = vv v + wv z ) 

[ V 2 - d t ]w - ^ = i( vw v + ww z ) 


(2.11a, b,c,d, e ) 


u = v = w = 0, y = 0,1. 

with V\ = d* + d\. 

The solution of the linearized form of the above equations shows that Taylor vortices 
grow exponentially in time for T > T e ~ 5162. and weakly nonlinear theory (Seminara 
(1975)) can then be used to show that the most unstable mode is stabilized by nonlinear 
effects. The neutral curve for the linear stability problem is shown in Figure 1. 

In the fully nonlinear regime possible stable finite amplitude solutions of (2.11) can be 
found by stepping the equations forward in time until the flow equilibrates. Bennett and 
Hall (1988) investigated the instability of these vortex flows to small Tollmien-Schlichting 
waves. Here in contrast we allow the Tollmien-Schlichting waves to be sufficiently large 
that they have an 0(1) effect on the vortex flows. 

Second, if the vortex wavelength is 0(J2e - ^) we define 

Z = Re~*z, 6 = -TRe~ 1 r, (2.12a, b) 

4 

and (2.10) is modified to give 

(2.13a, b) 

[u,v,«;] = [ tr o + “’ 2 ^’^rn i + 0 (' Re “^)]* 
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The functions u, v, w and p are now functions of y, Z, and t and satisfy 

dv dw 

!hj + dZ =0, 

i d l - d t}u = \vU' 0 + ]-{vu v + wu z ). 


1 - dp 

-TuU 0 - / = 0, 

2 dy 


(2.14a, b, c,d, e) 


{' d l ~ d t ~ H = \{™y + ww z ), 


u — v = w = 0, y = 0, 1. 

The linearized form of this system has unstable disturbances for any nonzero value of T. 
Indeed the neutral value of T is proportional to k~ 2 if k is the axial wavenumber. This 
means that finite amplitude solutions of (2.14) cannot be found by integrating forward in 
time since the energy of the disturbance will rapidly cascade into the higher harmonics. 
We return to the long wavelength limit in §4, since we concentrate next on the nonlinear 
interaction problem for = 0(1). 


3. The nonlinear evolution equations for TS waves and Taylor vortices of 0(1) 
wavelength. 

Since the only major difference between the present calculation and that of Bennett and 
Hall (1988) is that the TS wave now has an 0(1) effect on the vortex only the essential details 
of the flow structure will be given. The size of the TS wave is fixed by the requirement 
that the nonlinear terms in the y and z momentum equations be comparable with those 

in (2.7c,d). This forcing is most important away from the viscous wall layers of thickness 
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Re~’’ in which the TS waves adjust to the no-slip condition at the wall. It is convenient 
at this stage to define the small parameter e by 

c = Re~$. 

Suppose then that TS wave is proportional to 

E = exp[i{aex - f 

where the wavenumber a and the slowly varying frequency fl are real. The frequency 
and amplitude of the wave must vary on the vortex time scale in order to allow for the 
situation when the vortex flow is evolving in time. If the vortex flow is in equilibrium then 
the frequency and amplitude of the TS wave are constant. 
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Away from the viscous layers near y = 0, 1 the velocity field and pressure expand as 


f ( Uq + uo ^ 


eV ) 

E + complex conjugate 1 + 0(e 7 ) >. 

U={\ 6 7 /2 

+ 

6 7 V! 

V € 7 Wo/2 ) 


e 7 wi ) 

JL JJ (3.1a, b) 


p = (—12 x/Re + p 0 /2Re 2 + [e*piE + complex conjugate] + 0(e 7 ). 

The above expansions are substituted into (2.7) and the zeroth order system for the x 
independent part of the flow satisfies the vortex equations (2.11) with (u,v,tv, p) replaced 
by (uo, v 0 , wo, Po) and with forcing terms F\,F 2 on the right hand sides of (2.11c,d). These 
forcing terms come from the interaction of the TS wave with itself and are given by 


F\ = 2[(— iau\v\ + v\v* ly + tvivJJ + complex conjugate], 
F 2 = 2[(— iauiw\ + Viw[ y + wiw{ z ) + complex conjugate]. 


(3.2a, b) 


(In these two equations, * denotes complex conjugation.) The boundary conditions for 
(uo, Vo, wo,po) come from the no-slip condition 


u Q = vo = wq = Oony = Oand 1. 


together with a periodicity condition in z. This follows from the fact that the forcing terms 
become negligible in the wall layers. 

Now consider the zeroth order approximation to the TS wave equations in the core. 
We find that they are identical to the corresponding equations in Bennett and Hall’s linear 
analysis, so that writing U = Uq + U 0 we have 


taui + Vi v + w lz = 0, iaXJui + viU y + w\U z = 0 


iaU vi = —piy, iaU wi = — p lz , 


(3.3a, 6, c,d) 


together with the tangential slipping condition at the walls 


vi = 0, y = 0, 1. 


Bennett and Hall simplified these equations by introducing a normalized pressure term <j> : 

Pi = Pio{t) + ( ia) 2 A(t)<f>{y,z ) 

«i = AV 2 <f>/2U, vi = -iaA<t> v /U, w t = -iaA<f> t /U, 
where pio(t) and the amplitude A{t) need to be determined and <j> satisfies the core equation 

+ {jb), = °- ^Periodic inz. 

U U (3.4) 


<f> = 0ony = 0. <^ = lont/ = l. 
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A related coreflow problem arises in Smith (1977). We can now eliminate uijiqand Wi from 
the forcing terms in the vortex equations to give 


•(*.)* + (*.)* 


Fl = 2ot 2 \ A | 2 {[ ™ ]{ U\ + [ vm , : - ^ ] 


M 2 + (t*Y 


Fi = 2a 2 


t/ 4 

UvY + (^) 2 


u 2 


A J, ]{ u 2 ) z + 


r (^v ) 2 + (^g) 2 l 

1 c / 2 (<^) 2 J ‘ j 


(3.5a, 6) 


We note that the last terms in these equations can, if necessary, be absorbed into the pres- 
sure derivatives in the y, z vortex equations. In the two-dimensional case U = U(y),<f> = 
<j>(y) and we see that F 2 = 0 , and (using equation (3.5)) we note that Fi can be completely 
absorbed into p 0v . Thus for flows that are strictly two-dimensional nonlinearity occurs only 
at larger values of A, as in Hall and Smith (1982). 

Because the leading order wave equations are inviscid it is necessary to consider the 
flow within two viscous wall layers at y = 0, 1. The flow in these regions is exactly the 
same as in linear theory since the TS waves are sufficiently small. Thus 

y = e 2 Y 

in the lower viscous layer and 

u - {(c 2 A 0 ( 2 :, t)Y . , -e u no Z (z, t)Y 2 / 2, e 9 p 0 (z, t)T) 


(3.6a, 6) 


+\(e e U,e 9 V,e 7 W)E + complex conjugate ] H }, 

p* = pU 2 m {-12e 7 x-e™nozY 

+[(e 8 P + e 10 QE + complex conjugate ] + •••}. 

where A 0 can be determined from X 0 (z,t) = U v \ y =o. The equations to determine U,V,W,P 
are 

iaU + V v + W z = 0. 


*(-n + a\oY)U + A 0 V + A 0z YW = -iaP + U vv , 


(3.7) 


p v = q v = p, = o, t (-n a A 0 y)w = -q z + w vy , 

with boundary conditions U = V = W = 0 on Y = 0, together with matching condi- 
tions with the core : 


U — ► A(t)<py VV \y-o/{2X 0 ),W — ► 0 and P -* p w as Y -*■ oo, 

and a periodicity condition. An important feature of these equations is the presence of a 
second order pressure term Q. We need to include this term because the z variation of the 
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flow is relatively fast compared to the x variation of the waves, 0(1) as opposed to 0(c), so 
the 2 -derivatives of this smaller pressure will come in at higher order. 

These boundary-layer equations were solved by Bennett and Hall (1988) in terms of 
Airy functions Ai(f) to give an eigen- or dispersion-relation. The equations governing the 
wall layer flow are then matched to the coreflow solution to give the following first order 
differential equation in Q z , 

Qz. t + = <* 2 Pio - \<x 2 Mv VV \v=o — ( 3 * 8a ) 

Ao * (tG£Ao)SA(£) 

where 

m = -(§ + + A.-'(f)]), (3.84) 

QpW 6 roo 

£ = j- andif(^) = / Ai(s)ds. (3.8c, d) 

(aA 0 )* J i 

A similar equation involving Ai( 2 ,t) = U v \ v =\ arises from the boundary-layer at y = 1. By 
applying the periodicity condition in z, combining the two results to eliminate p 10 and A 
we can determine the eigenrelation. This is somewhat unweildy to write down but can 
best be calculated from the following: 


2 * y 


(3.9a) 


where 2 tt/A: is the period of the vortices and f n { z ) and g n { z ) satisfy the following first order 
differential equations, 


f,. + 

A n (*o:A n )3ff(^ n ) 


9nz + = 1, 


with boundary conditions 


/n(o) = 0, flf n (0) = 0, 

with V’ and K defined above and £„ given by 

He *'/ 6 

' . 3 , « = 0,1. 

(aA„)s 


(3.96, c) 


(3.9d) 


(3.9e) 


In order to facilitate the numerical solution of the vortex equations it is convenient to 
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eliminate the pressure from the y and z momentum to give 

{V 2 - d t }u 0 - \vqUq = Mu 0 , 

{V 2 - d t }V 2 v 0 + TU 0 u 0zz = (Mv 0 ) zz - ( Mw 0 ) zz - iT(u 2 ) s 
+2a 2 1 A | 2 


(3.10) 




«0 = V 0 = WO = 0, y — 0, 1, 

where M = \{vod v + w Q d z } and 0 appearing on the right hand side of the equation for v 0 
satisfies the ’core’ equation (3.4). 

In the absence of a Tollmien-Schlichting wave it is well-known that, without any loss 
of generality, we can take u 0 , v 0 to be even functions of z and w 0 is then an odd function 
of z. In fact (3.4), (3.10) and the wall equations show that the forced equation also admits 
a solution with tt 0 5 v o even in z with w 0 odd in z. In this case Q and the corresponding 
upper wall pressure together with <f> are even functions of z. Thus at 0(1) wavelengths the 
mixed vortex-TS state is found by solving (3.10) subject to the TS wall equations (3.9b, c). 


The numerical method. 

When the Tollmien-Schlichting waves are small, as in the linear problem, the terms 
Fi and F 2 become zero to leading order and the vortex and dispersion equations then 
decouple. Bennett and Hall solved equation (3.10) with Fi = F 2 = 0 using the method 
of Rogers and Beard (1969) for solving the Taylor vortex equations, by expanding the 
vortex velocities in Fourier sine and cosine series in z and by advancing in time to reach a 
steady solution. The nonlinear terms were calculated explicitly. Equation (3.4) was then 
solved for the wave pressure <f> by using finite differences in both y and z and iterating to 
a solution. Several thousand iterations were required, the exact number depending on the 
method used and the step lengths in y and z. Once the skin-frictions A 0 , Xi and pressure 
derivatives 4> yvv \ y =o t \ had been determined a and fl were found by solving the dispersion 
relation (3.9) using Runge-Kutta methods. The Airy functions were calculated using a 
combination of ordinary and asymptotic series, depending on the value of the argument £. 

We could adapt this approach to solve our unsteady nonlinear problem in the following 
way. First we fix the T-S wave number a, the initial vortex velocity Uo(to) and wave 
pressure <£(t 0 ). We then guess an amplitude at the next time step | A{ti) | and use it and 
<f>{to) to calculate the forcing terms F\ and F 2 . The vortex equations (3.10) can then be 
solved for «o(ti) while <f>(ti) can be found from the pressure equation (3.4) using the new 
velocity. The eigenrelation (3.9) can then be solved for H(fi). The guess A(fj) can then 
be iterated on, until the imaginary part of fl is zero. This process could then be repeated 
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at the next time step t 2 . We would then be able to observe the wave amplitude | A(t) \ 
evolve in time for various wave-numbers a away from the neutral point. The disadvantage 
with this method is that at each time step equation (3.4) needs to be solved accurately 
several times. The iterative method used by Bennett and Hall to solve (3.4) was very slow, 
taking up about half of the computing time for the entire linear calculation. We therefore 
feel a quicker scheme is needed for this part of the problem in order to carry out extensive 
numerical investigations using this approach. 

It is possible however to calculate solutions to the steady nonlinear interaction problem 
less expensively. The method used is as follows. First we fix the term a 2 | A | 2 that occurs 
in equations (3.10). Starting with a guess for and <f> we march forward in time with the 
uo, Vo equations. At each time step we update <f) from equation (3.4) by performing a few 
iterations. The forcing terms Fi and F 2 can then be updated for the next time step. We 
continue this process until u$ and <f> have settled down to steady solutions. The dispersion 
relation (3.9) can then be solved for real a and Cl. Once a is known the steady amplitude 
| A | can be determined. Even if we were to solve the pressure equation (3.4) accurately 
at each time step, the velocities «o(t) that are produced before a steady solution emerges 
are not solutions of the physical problem. This is because we would find when we solved 
the eigenrelation (3.9) that in order to make Cl real we would have to have different a’s at 
each time step, whereas the theory requires a to be independent of time. Only when the 
calculation has settled down to a steady solution do we have anything that makes sense 
physically. This means that the algorithm may well converge on unstable solutions as well 
as the stable ones, although of course we will not be able to tell which is which. Because <f> 
is only iterated upon a few times at each time step the overall time taken for this method 
is about equal to that for the linear calculation. 

Since the TS wave part of the equations is still linear the calculation of <f> and the 
solution of the eigenrelation (3.9) are unchanged from Bennett and Hall. The algorithm 
used to step forward (3.4), (3.10) in time is virtually identical to that used by Bennett and 
Hall. The only difference is the presence in (3.10) of the forcing terms F x and F 2 which 
are periodic functions of z respectively. All the periodic functions of z can be expanded in 
the form ^ 

««=i: ««.«*•** ( 3 - n ) 

-oo 

where k is the vortex wavenumber. We recall that the equations from the interaction 
allow a solution with uo, t>o and the TS wall pressure even in z and wo, po odd in z. Thus 
we expect that there will be a solution of the nonlinear problem with uo n = uo-n> etc. 
However, in order to allow for the possibility of solutions without this symmetry we did 
not force this constraint. A consequence of this is that there is a mean (independent of 
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x, z) flow induced in the axial direction and an induced pressure gradient in that direction 
is needed to satisfy a mass flow constraint. A similar pressure function arises in the work 
of say DiPrima and Stuart (1975) or Hall (1984). Thus we allowed for such a function in 
our calculations. However, the only nonlinear solutions found had the above symmetries in 
z in which case there is no mean axial flow. It is possible that nonlinear solutions without 
these symmetries exist but they apparently do not bifurcate from the pure vortex flow. 
We postpone until §6 a discussion of the numerical results obtained following the scheme 
outlined above. 


4. The strongly nonlinear interaction between long wavelength Taylor vortices 
and Tollmien-Schlichting waves. 

Here it is assumed that the channel curvature parameter 6 defined by (2.1) is sufficiently 
large to support Taylor vortices with cross-stream wavenumbers of order e = Re~% More 
precisely we consider the limits 6 — ► 0,Re — ► oo with T defined by (2.12b) held fixed. In 
this case the basic circumferential flow can support vortices and oblique TS waves with 
comparable axial wavenumbers. The TS waves now have the structure first discussed 
by Smith (1978) in connection with the instability of the flow in an elliptical pipe. The 
overall size of the TS wave is again fixed by the condition that it should be sufficiently 
large in the core to drive the vortex flow at zeroth order and indeed alter the mean-flow 
profile. The streamwise and time dependence of the wave will again by expressed in terms 
of E = exp i{aex — f* dt } in order to account for the possible evolution of the TS 

frequency as the vortex flow develops in time. The axial dependence of both modes is 
now of course entirely on the Z = tz scale. In this case the appropriate expansions of the 
velocity and pressure fields now become 



f U o + uo ^ 


e 5 u x \ 


e 7 v 0 /2 

+ ••• + 

e 6 vi E + complex conjugate H ► 


i e 6 «>o / 2 ) 


< e 7 «h J 


(4.1a, 6) 

p = -12xe 7 + p 0 Sj- H 1- JpiE + complex conjugate -I . 

Here again the Taylor vortex functions Uo, t>o, iy 0> Po and the TS functions «i, t>i, tui, pi must 
be matched with the appropriate expansions in the wall layers. In fact the Taylor vortex 
functions satisfy (2.14) but with forcing terms arising from the interaction of the Tollmien- 
Schlichting wave with itself. The zeroth order core equations for the TS wave yield the 
solution 

ui = AU V , Vi = —iaAU, pi = — Aa 2 f? U 2 dy, 

2 

(4.2 a,b,c,d) 

gkgg * l 


W 1 = 


iaU 


12 



Here A is an amplitude function of Z = ez and the Taylor vortex time scale t. It should be 
noted that A cannot vary on a faster time scale more appropriate to a TS wave alone since 
the Taylor vortex cannot itself respond on such a scale induced by the forcing terms. With 
the form of the TS wave determined in the core the forcing terms in the vortex equations 
in the core can be expressed in terms of A and U. If the vortex pressure is eliminated we 
obtain 

<9t>o | drop rj 

dy dZ ~ 

{dy - dt}u 0 = \v q Uq v + Mu 0 , 

{dl " dt}d}v o + TUoUozz = -{Mw 0 ) zz - \T{ul) zz + 2**{\ A || ( U 2 ) v } z , 

(4.3) 

u 0 = v 0 = w 0 = 0, y = 0, 1 


where 


X = v 0 d v + w 0 d z . 

We have anticipated above that the vortex flow in the presence of a TS wave should still 
satisfy the no-slip condition at the wall. In fact the axial momentum equation for the 
vortex flow in the wall layers is forced by the TS waves. However the forcing function 
decays to zero like the inverse square of the distance from the wall so that the appropriate 
matching conditions for wo are that it should vanish at y = 0, 1. 

It remains for us to determine the structure of the TS wave and Taylor vortex in the 
viscous wall layers. The appropriate expansions in the wall layer at y = 0 are 

u= {[ £ J Ao(Z ) <)y, ( ,1 K 0) £X] 

[(e 5 U,e s V,e s W)E + complex conjugate] + ■ • •} (4.4a, b) 

p = — 12c 7 x + e 12 P 0 H b ( e 1 PE + complex conjugate ) H 

The axial momentum equation for the vortex is now the only equation forced by the TS 
wave. However, the equations for the TS wave depend only on the vortex flow through 
the shear A 0 so we do not solve for V 0 ,W 0 and P 0 here. It suffices to say that the vortex 
equations in the wall layer at y = 0 can be solved such that the solution matches with the 
coreflow vortex solution. The zeroth order equations for the TS waves in the wall layer at 
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y = 0 are 


ia.U + Vy + Wz 


= 0 


i(-fl + aX 0 Y)U + A 0 V + A 0 Z YW = -iaP + U YY 


Py 


= 0 


—Pz + W YY 


i{-n + aX 0 Y)W 
with boundary conditions 

U = V = W = 0, Y = 0 

and matching condition 

U — ► A 0 A, Y —* oo. 

Following Smith (1979) we can solve (4.5d) to give 


W = 


Pirn 


* 2 5 

As 


where 


£ - A *tt) f( 0 7$?} ^.(5i)ds, 

The x velocity component of the TS wave is then given by 

A*£7 e XoPzzC'U) + lA„ z P z {3£'({) + j£“(f) 

which matches with the coreflow solution if 

AlA„X 1 £’(&) 

Jf(fo) 2 ozz “if({„)’ 


(4.5a, 6, c,d) 


(4.5e) 

(4.5/) 

(4.6) 


(4.7) 


(4.8) 


with K(£o) defined by (3.8d). A similar expansion procedure can be carried out in the 
upper layer, while the core equations for the TS wave show that P and Q , the pressures 
in the lower and upper layers respectively, are related by 


P — Q = a 2 A f 1 U 2 dy. 

Jo 

The TS pressure P can be written using (4.5b) as 


(4.9) 


iaP = A 3 U(£ 


l«=fo 
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so that (4.7), (4.8) now give 


» _i_ ^0./.. D> _ -^<(^o)Ao[Q + a 2 p 


P" + ^-tpoP' = 
Ao 


K 0 ct 2 J 


where tp 0 is defined by (3.8b) with £ = £ 0 and 

J= f 1 U 2 dy. 

Jo 

A similar calculation in the upper layer shows that 


(4.10) 


(4.11) 


Q" + %hQ' = o?Q + Q] . (4.12) 

Ai KoorJ 

Here Ai = — | u= i and £i,V>i are obtained from £ 0 , 4>o by replacing A 0 by Ai in the 
definition of these quantities. Thus the evolution equations (4.3), (4.9), (4.10), (4.12) 
determine the TS and vortex structure as the nonlinear interactions take place. The 
vortex velocity components uo, Vo, w 0 , the Tollmien-Schlichting amplitude A (or P, Q) and 
frequency fi are all functions of the vortex time scale t. Again it should be noted that 
the vortex structure cannot tolerate an evolution on a faster time scale typical of TS wave 
growth. 

Suppose then that, for a given ot,T, we are given u 0 , v 0 , w 0 , P, and fl. If a guess is 
made for | A \ max at the next time step (4.3) can be advanced using some finite difference 
discretization in t and y together with a Fourier expansion in Z. The wall equations (4.10), 
(4.12) at the new time step constitute an eigenrelation fl = f!(a) and this can be solved for 
the complex fl which satisfies this eigenrelation. It is however necessary in our expansion 
procedure for fl to be purely real so the calculation is repeated with a different | A \ maz in 
order to make fl,- = 0. At each stage of this iteration procedure A(z) is written as 


A — A max Jf ( 2 ) 


and A(z) is taken from the solution of the eigenrelation at the new time step. When 
the iteration has converged the scheme can be used to advance the solution again. If 
there are stable equilibrium states ( d t = 0) then we would expect that after a sufficiently 
long time the vortex and TS wave would become steady. The above procedure will then 
follow the time dependent evolution of the vortex-TS interaction. Because this calculation 
must be carried out over a significant time interval and there are many values of a,T 
to be investigated it was decided to concentrate on an alternative scheme which would 
locate any equilibrium states directly. In this scheme d t is set equal to zero and values for 
a, | A \ max , T are chosen and a guess is made for A(z). We then use a Newton iteration 
procedure on the steady vortex equations to find the corresponding equilibrium vortex flow. 
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The TS wall equations are then solved for the real a and fl which satisfy the eigenrelation. 
The new amplitude function A (z) obtained from this calculation is then used in the vortex 
equations to find an ’updated’ vortex flow. The wall equations are then used to ’update’ 
a, Cl and the procedure continued until the vortex solution and (a, fi) converge for the 
given value of a, \ A \ max and T. 

It was found that the above iteration procedure converged and typically about 100 
steps in (0, 1) were taken to integrate the vortex equations using a Runge-Kutta scheme. 
The wall equations were solved by expanding P,Q and the coefficients in (4.10), (4.12) 
in Fourier series in Z. It was found that about 8 Fourier terms were usually adequate to 
solve the eigenrelation to the accuracy needed to plot the figures given later in this paper. 

5. Secondary instability of two-dimensional Tollmien-Schlichting waves de- 
scribed by the TS-vortex interaction equations. 

Here we demonstrate how the interaction equations can be used to model the three- 
dimensional breakdown of two-dimensional TS waves in a straight channel. Suppose then 
that A(Z,t), the wave amplitude in the core region, takes the form 

A(Z, t) = Ao(t) + Aie tut cos /3Z (5.1) 

where Ai is independent of time and is small compared with the 2D amplitude Aq. If 
A\ = 0 the TS forcing term in (4.3) vanishes identically so that Aq does not itself induce a 
longitudinal vortex field. This means that Ao(t) evolves as the solution of a cubic amplitude 
equation which must be found at higher order. In fact A 0 (t) evolves on a longer time scale 
than t so that in (5.1) Ao can be treated as a constant. We shall determine the values 
of Ao at which the 3D perturbation in (5.1) is neutrally stable. The precise form of the 
perturbation in (5.1) corresponds to two equally inclined oblique TS waves of the same 
amplitude. 

The vortex velocity field (uo,vo,wo) induced by the TS wave (5.1) can be expressed as 
(u 0 ,v 0 ,w 0 ) = [4a 2 fl 2 A 1 A 0 (u 0 i,voi,woi) cos j3Ze ,lJt + complex conjugate] + 0(Aj). (5.2) 

Here (uoi,uoi>woi) satisfy the problem 


# v oi a 

a 01 
oy 

= 0 


(dj - *w)u 0 i 

= 

(5.3) 

(d 2 y - iu)d 2 y v 01 

= - Wo 


UQl = Voi = v oi = 0, 

y = 0,1 
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which must be solved numerically (see below) . 

In the lower and upper wall layers P and Q expand as 

P = P 0 + P 1 cos/3Ze i “ t 


Q = Qo + Qi cosfiZe 


tUJt 


and from (4.9) it follows that 


Pi~Qi = + SAiA 2 0 a 2 p 2 Ji] 


where 


J\ = [ U 0 u 0 idy. 

Jo 

In fact it can be shown that Qi = —Pi. If the wall equations are expanded in powers of 
A, then at zeroth order we obtain 


a 


(5.4) 


2 _ 5.6f>i; 0 

K 0 A? 

which yields the usual neutral values of a, Cl when 

— ~ 1.001, £ = — 2.298t*. 

X 

The condition (5.4) ensures that the two-dimensional TS wave is in neutral equilibrium. 
At next order we find that 

20 

T 

where 


1 I — tf 2_2 a 2 t 2 a 2 r ® r , U Ol( 0 ) 

1 + —p a A 0 J 1 — a A 0 (— Ji -I — 


*] 


(5.5) 


_5_ _ ZAjp , £oi J_ 

18 9 K A i0 K 0 h 

The integral Jj and the shear «oi(0) depend on w and (5.5) can be solved by taking real 
and imaginary parts to give 

/? 2 = 86/d, 
aM 2 - * 

0_ «Ffi 


where 


(5.6a, 6, c, d) 


5 $uL(0) 

“6 Jl+ 2 _2 = “ + ‘ 6 ’ 
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Ji — c + id. 


17 


A numerical investigation of (5.3) shows that (5.6a,b) give / * 2 > o, ol 2 AI > 0 for ui > 59.5. 
In fact the right hand side of (5.6b) is positive for u> > 0 while the right hand side of (5.6a) 
— * oo when /?a -1 — > oo, and so 3D waves oriented at all angles to the flow direction can be 
neutral. However, the value of A 0 corresponding to /? increases monotonically with /?a -1 
so that 3D waves propagating almost parallel to the flow direction (i.e. with small (3 ) are, 
in a sense, the most dangerous. Thus (5.6a,b) determine the amplitude of the 2D TS wave 
which is neutrally stable to a pair of oblique waves with wavenumber /?. It follows that at 
this amplitude there can be a secondary instability of the 2D TS wave to a pair of equally 
inclined oblique modes. 

6. Results and Discussion. 

We first discuss our results for the interaction problem of §3 which we recall concerns 
Gortler vortices of wavelength comparable with the channel width. In this case, finite 
amplitude vortex flows are possible only for T > T c = 5162 and the least stable ax- 
ial wavenumber has k = 3.951. Since there are no unstable TS waves with streamwise 
wavenumbers 0(Re~%) and spanwise wavenumbers 0(1), we restricted our calculations to 
the case T > T e and spanwise wavenumber k = 3.951. 

In Figures 1, 2, we show the dependence of a \ A | and fl on a for equilibrium states 
corresponding to T = 11,000, and T = 22,000 respectively. We see that at these Taylor 
numbers the neutral wavenumber and frequency both decrease monotonically as | <xA | 
increases. Results for T = 5500 have the same properties and are shown in Table 1. From 
these results, we deduce that the bifurcation diagram for a fixed frequency disturbance is 
as sketched in Figure 3. We have not computed sufficient results to draw such a figure 
in detail because it would require an order-of-magnitude more computing time than that 
needed to generate Figures 1, 2. Linear theory shows that for each of the points of Figures 
1, 2 the flow with infinitesimal | A | is unstable. We conclude that the interaction problem 
at 0(1) vortex wavelengths leads to a supercritical bifurcation to a stable mixed TS-vortex 
state. 

Now let us turn to our results for the interaction problem specified by (4.3), (4.9), 
(4.10) and (4.11) corresponding to a longer-scale spanwise dependence. In the absence of a 
finite amplitude Gortler vortex, we find that the neutral values of a and fl corresponding 
to an oblique TS wave with A ~ cos z or A ~ sin z are 

a = 4.396 , n = 20.456. 

In the first instance we confine our remarks to the situation when T > 23711 so that a 
finite amplitude vortex can exist in the absence of TS waves. In Figure 4 we show how 
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a, ft vary in the presence of a finite amplitude vortex at different values of the Taylor 
number. These values correspond to the linearized problem A — *■ 0 and we see that both 
a and ft increase monotonically with T. We note that a constant frequency TS wave will 
change from being stable to unstable when T is decreased through any point on the curve 
of Figure 4 for T greater than its linear neutral value ~ 23711. 

Next, suppose that the Taylor number is held fixed and the neutral values of a and ft 
are calculated for different values of | A \ max . When | A | mai — ► 0, then the neutral values of 
a and ft must tend back to the neutral values appropriate to the current value of the Taylor 
number. The results of such a calculation are shown in Figure 5 at two typical values of 
the Taylor number. We see that a. and ft increase monotonically when | A \ max increases. 
We conclude that the bifurcation picture for any constant frequency TS wave will be as 
sketched in Figure 6. At first sight, it appears that this is a subcritical bifurcation to a 
mixed vortex-TS state. However, we found earlier that a constant frequency linear TS 
wave changes from being stable to unstable when T is increased, so that, if we associate 
the terms ’subcritical’ and ’supercritical’ with the linearly stable and unstable regions, we 
conclude that then a constant frequency disturbance undergoes a supercritical bifurcation 
to a finite amplitude state. A weakly nonlinear analysis of (4.3), (4.9), (4.10), (4.11) shows 
that when this bifurcation occurs the finite amplitude solution is stable and we do not 
expect that there will be any finite time singularities of the full time-dependent problem 
associated with the finite amplitude state. 

Now let us consider the situation when T < 23711 so that in the absence of a Tollmien- 
Schlichting wave there is no vortex activity. In order to see the differences which emerge in 
this case, we shall now indicate briefly how an amplitude equation typical of those obtained 
for weakly nonlinear stability problems using the Stuart- Watson method can be retrieved. 
Suppose then that the Taylor number T is held fixed and that the neutral values of a and 
ft are given by 

ot = (Xff , ft = ftw* (6.1) 

These are the neutral values appropriate to | A | mai — ► 0 , so that if we now write 

| A \ max = 6 J 


where S is small and positive we anticipate that with T held fixed the appropriate expan- 
sions of a and ft become 

a = ajv + <5 a + • • • , 


ft = ftjv + 5ft + • • • . 


(6.2a, b) 
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We see from (4.3) that u 0 , v 0 , w 0 in (4.3) must be expanded as 

u 0 = 8u cos2z H , 


where u,v,w satisfy 


v 0 = 8v cos 22 + • • • , 
w 0 = 8w sin2z + • • • , 


d?u 

dy 2 


W V ^ ^ ^ 

— + 2w = 0, 
dy 



d 4 v 

dy 4 


4TU 0 =-4 a 2 TUS y , 


(6.3a, b) 


u — v = tv = 0 , y = 0, 1. 

We note that the homogeneous form of (6.3) has a solution if 4 T ~ 23711, so that 
u,v,iv then become singular near T = T + ~ 5978. In fact, these functions behave like 
j r zf+ so ^at u,v,w necessarily change sign at T + . The wall equations (4.10) and (4.12) 
together with (4.11) are then found to have a solution only if a,/3 satisfy an equation of 
the form 

a[c + id] + f2[e + if] + g + ih = 0. (6-4) 

Here c,d,e,f,g and h are real constants , the first four of which arise from the expansions 
of a and fi in (4.10) - (4.12) , whilst g and h arise from u , the order 8 correction to 
the coreflow. By taking real and imaginary parts in (6.4) we can solve for £1 and a. Our 
numerical solution of the full problem for T > 23711 suggests that there 0 is always 
positive. However, because of the singularity of u at T + it follows that 0 must change sign 
at f = f + . 

It follows from the above that in the neighbourhood of T + there is a change in the 
bifurcation structure of the mixed vortex-TS state; thus on one side of this Taylor number 
the bifurcation will be supercritical whilst on the other side it will be subcritical. In 
fact, this structure is discussed by Hall and Smith (1988a). In the present investigation 
the iteration scheme used to solve the interaction problem of §4 failed to converge for 
Taylor numbers less than 25000. Until the method failed to converge it was found that the 
bifurcation to the mixed state was always supercritical. This suggests that at sufficiently 
high Taylor numbers the time dependent interaction problem does not have a finite time 
singularity but that at lower Taylor numbers there could be such a singularity associated 
with the subcritical bifurcation. 
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Fig. lb. The dependence of |ai4| on fi for T = 11,000, k = 3.951. 
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Fig. 2a. The depend 
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Fig. 4. The dependence of a and n on T. 
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Fig. 5a. The dependence of |.A| max on a for the case T = 28511. Fig. 5b. The dependence of |yi| mai on fi for the case T = 28511. 
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Fig. 5c. The dependence of on a for the case T = 31000. Fig. 5d. The dependence of |v4| mai on fi for the case T = 31000. 
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Fig. 6. A sketch of the bifurcation diagram for a constant frequency TS wave. 
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